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Summary 

This document describes the initial formulation 
(Version 1.0.0) of the Atmospheric Radiation Measure- 
ment (ARM) Program satellite data analysis procedures. 
Techniques are presented for calibrating geostationary 
satellite data with Sun synchronous satellite radiances 
and for converting narrowband radiances to top-of-the- 
atmosphere fluxes and albedos. A methodology is docu- 
mented for combining geostationary visible and infrared 
radiances with surface-based temperature observations to 
derive cloud amount, optical depth, height, thickness, 
temperature, and albedo. The analysis is limited to two 
grids centered over the ARM Southern Great Plains Cen- 
tral Facility in northcentral Oklahoma. Daytime data 
taken April 5 to May 1, 1994 were analyzed on the 0.3° 
and 0.5° latitude-longitude grids that cover areas of 
0.9° X 0.9° and 10° x 14°, respectively. A variety of 
cloud conditions ranging from scattered low cumulus to 
thin cirrus and thick cumulonimbus occurred during the 
study period. Detailed comparisons with hourly surface 
observations indicate that the mean cloudiness is within a 
few percent of the surface-derived sky cover. Formats of 
the results are also provided. The data can be accessed on 
the Internet via the World Wide Web at the following 
uniform resource locator: 

http: //albedo. larc . nasa.gov: 1123/arm.html 

Introduction 

The Atmospheric Radiation Measurement (ARM) 
Program is a long-term measurement and modeling pro- 
ject designed to improve the understanding of atmo- 
spheric radiation and its interaction with the atmosphere 
and surface (Stokes and Schwartz 1994). The program 
concentrates on determining the influence of clouds and 
their radiative feedback effects on climate. The primary 
objectives are (1) to relate observed radiative fluxes in 
the atmosphere to the atmospheric temperature, composi- 
tion (particularly water vapor and clouds), and surface 
radiative properties and (2) to develop and test para- 
meterizations of atmospheric water (clouds and vapor) 
and the surface characteristics affecting atmospheric 
radiation. The parameterizations are intended for use in 
prognostic mesoscale and general circulation models 
(GCM). The models and measurements produced by 
ARM will be extremely valuable resources for advancing 
our understanding of the role of clouds in climate 
change. 

The ARM observational program focuses on contin- 
uous long-term measurements taken from a variety of 
sensors at several locales representing various climatic 
regimes. More specifically, a complete suite of surface 
instruments is situated at a central facility in a given 
locale. Other sites surround the central facility at dis- 


tances sufficient to monitor portions of an area equiva- 
lent to a GCM grid box. These extended facilities operate 
with a reduced complement of sensors. The locale and 
the regions around it are also monitored by Sun synchro- 
nous and geostationary meteorological satellites. The 
operational surface and satellite measurement system is 
occasionally enhanced by aircraft and additional surface- 
based instrumentation during special intensive observing 
periods (IOP). The first ARM site located in northcentral 
Oklahoma, called the Southern Great Plains (SGP) 
locale, will be followed by other sites in areas such as the 
tropical western Pacific and Alaska. The long-term ARM 
measurement systems are evolving at each locale as new 
instruments are deployed and the latest analysis tech- 
niques are implemented. The ARM surface-based sen- 
sors and the products derived from these techniques have 
been described in several reports (e.g., Schneider, Lamb, 
and Sisterson 1993). 

The ARM satellite measurements complement the 
surface and atmosphere observations. Measurements of 
radiative fluxes at the top of the atmosphere and the sur- 
face can be used to determine the energy budget of the 
atmosphere. Many ARM surface instruments provide 
continuous, high-resolution measurements of various 
properties over a small area within the extended locale. 
The satellite can yield contiguous, low-resolution data 
over the entire locale and its surroundings. Ideally, the 
satellite data can be used to extend the information 
gained by the surface instruments by relating the surface- 
observed quantities to those deduced from the satellite 
radiances. The resulting relationships can be applied at 
other locations observed only by the satellite to infer 
quantities that are normally measured from surface 
instruments. To be useful for such applications, the radi- 
ances measured by the satellites must be analyzed to 
derive values for parameters that are relevant to cloud 
and radiation processes. 

This report documents the initial analysis procedures 
that will be applied to the ARM satellite datasets; these 
procedures are collectively designated Version 1.0.0. 
Examples of the initial products are also included with 
preliminary validation. The ARM satellite data analysis 
program will ultimately be applied on an operational 
basis to geostationary and Sun synchronous datasets at 
all hours of the day. To reach this objective, the coverage 
and complexity of the analysis algorithms will be incre- 
mentally increased. Daytime analyses of geostationary 
satellite data for limited periods will be followed by con- 
tinuous daytime analysis. The first digit in the Version 
1.0,0 designator refers to the use of a visible-infrared 
bispectral method applied only during daytime. This first 
digit will denote changes such as the inclusion of new 
channels or nighttime data. The second digit in the desig- 
nator will denote changes in the bispectral algorithm, and 



the final digit will indicate changes in the calibrations 
applied to the satellite data. Additional procedures for 
analyzing daylight Sun synchronous data will be imple- 
mented in the firture to include cloud effective particle 
size. Nighttime datasets will be analyzed on a periodic 
basis until a complete package of both daytime and 
nighttime procedures is obtained. The first step in this 
process is described here together with a summary of the 
results derived specifically for the Spring 1994 ARM 
IOP conducted over the SGP locale. 

Nomenclature 

ARM Atmospheric Radiation Measurement 

Program 

AU astronomical unit 

AVHRR Advanced Very High Resolution 

Radiometer 

a, b, c regression coefficients 

B Planck function 

C satellite-derived cloud fraction 

^ 'sfe surface-observed cloud cover 

CLD, CLR cloudy and clear, respectively 
CS cirrostratus model cloud 

D digital brightness counts, 0 to 255 

d s [ days since launch 

ERBE Earth Radiation Budget Experiment 

ERBS Earth Radiation Budget Satellite 

f D forward scattering truncation factor 

GCM general circulation model 

GOES Geostationary Operational Environmental 

Satellite 

g visible sensor gain 

h local time, hr 

HBTM hybrid bispectral threshold method 

IOP intensive observing period 

IR infrared, 10.5 to 12.5 pm 

ISCCP International Satellite Cloud Climatology 

Project 

L radiance, W/m 2 -sr 

LBTM layer bispectral threshold method 

LT local time 

LW longwave, 5.0 to 50.0 pm 

M ir , M lw infrared and longwave fluxes, W-nT 2 

McDDAS Man-computer Interactive Data Analysis 

System 


p pressure, hPa 

P m , Pfo P p temperature-reflectance pairs correspond- 
ing to clouds at 2 km, 6 km, and tropo- 
pause, respectively 

R linear correlation coefficient 

RH relative humidity, percent 

SCF SGP Central Facility 

SE standard error of estimate, percent 

SGP Southern Great Plains 

SW shortwave, 0.2 to 5.0 pm 

t a i transmittance of ozone above cloud 

t c i, t c T downward and upward cloud transmittance, 

respectively 

T equivalent blackbody temperature 

T a , T c , T t air, cloud-center, and cloud-top tempera- 
ture, respectively, K 

T cs , T s dear-sky and surface shelter temperature, 
respectively, K 

T lim , Tn m \, clear-sky temperature limits, K 

T liml 

T m , T h , T p temperature at 2 km, 6 km, and tropopause, 
respectively, K 

TAFB Tinker Air Force Base 

TO A top of atmosphere 

UTC Universal Coordinated Time, hr x 100 

VAFB Vance Air Force Base 

VIS visible, 0.55 to 0.75 pm 

VISSR Visible Infrared Spin Scan Radiometer 

WAFB Whiteman Air Force Base 

z altitude, km 

z c , z s , z t cloud-center, surface, and cloud-top 

altitude, respectively, km 

a albedo 

a c , a cd cloud and diffuse cloud albedo, 
respectively 

a cs , a sd dear-sky and diffuse clear albedo, 

respectively 

albedo of Rayleigh layer above cloud 
a w , otj, shortwave and visible albedo, respectively 
a swce effective shortwave cloud albedo 

y IR limb-darkening function 
8 Earth-Sun distance correction factor 

AD clear-sky count tolerance 
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A T scs greatest expected difference between T cs 

and T s , K 

A z cloud thickness, km 

e, z t cloud-center and cloud-top effective emis- 

sivities, respectively 

C, reflectance parameterization normalization 

factor 

0 viewing zenith angle 

0 O solar zenith angle 

p. cosine viewing zenith angle 

p 0 cosine solar zenith angle 

p reflectance 

p cl combined reflectance of cloud and over- 

lying Rayleigh layer 

p R2 direct reflectance of Rayleigh layer below 

cloud 

a standard deviation 

x cloud visible optical depth 

% 0 , ozone and Rayleigh optical depth above 

3 cloud, respectively 

X anisotropic reflectance correction factor 

<j> relative azimuth angle 

Subscripts: 

A AVHRR 

c clouds 

cs clear sky 

G GOES 

ir infrared 

Data 

The ARM Spring 1994 IOP operated for 21 days 
during April at the SGP Central Facility (SCF). Satellite 
data are analyzed for a 27-day period from April 5 to 
May 1, 1994. Other satellite data taken during April 1985 
are used to develop relationships for computing broad- 
band fluxes from narrowband radiances for the IOP 
dataset. Two grids are used for the ARM analyses. A 0.5° 
latitude-longitude grid, designated the mesoseale grid, 
extending from 32°N to 42°N and from 91°W to 105° W 
was selected to center the domain on the SGP, to include 
several GCM-scale grid boxes, and to minimize the diffi- 
culties involved in analyzing data over mountain snow 
fields. A much smaller fine-scale domain, extending 
from 36.16°N to 37.06°N and from 97.04°W to 97.94°W 
is divided into nine 0.3° grid boxes. This small-scale 
domain is centered on the SCF at 36.61 °N, 97.49°W 
(fig. 1). The fine-scale grid is designed for direct inter- 


comparison between satellite and surface observations. 
The 0.3° grid box is the minimum size needed for accu- 
rate cloud property retrievals. The mean surface eleva- 
tion z s was computed in kilometers above mean sea level 
for each grid box with the surface elevations given in the 
10' resolution U.S. Navy Surface Elevation Map (avail- 
able from the National Center for Atmospheric Research 
in Boulder, Colorado). The average elevations (fig. 1) 
range from a few hundred meters in the southeast to 
nearly 2.4 km on the western edge of the domain. The 
SCF is located at an elevation of =0.3 km. 

GOES-6 

Regional two-dimensional histograms were formed 
from April 1985, hourly 8-km visible (VIS) 8-bit digital 
brightness counts £>g and infrared (IR) equivalent black- 
body temperatures. The digital brightness counts and 
blackbody temperatures were measured by the GOES-6 
(Geostationary Operational Environmental Satellite) lo- 
cated over the equator at =108°W. Each region corre- 
sponds to one of the 2.5° latitude-longitude grid boxes 
between 32.5°N and 42.5°N and from 95°W to 105°W. 
Minnis, Harrison, and Young (1991) computed hourly 
cloud amount for each region from the individual radi- 
ances with the hybrid bispectral threshold method 
(HBTM). (See Minnis, Harrison, and Gibson 1987). The 
mean equivalent blackbody temperature T was computed 
from the IR radiances for each region and hour only for 
April 1 to 20 to eliminate a period of large calibration 
uncertainty that occurred during the last third of the 
month (Minnis, Harrison, and Young 1991). The narrow- 
band reflectance pg was computed with the calibration 
coefficients reported by Whitlock, LeCroy, and Wheeler 
(1994). Visible data were used for the entire month. 

For April 1985, the GOES-6 narrowband reflectance 
is as follows: 

p 6 = (0.0085Dg - 8.0) /526.2p o 5 (1) 

where \i Q = cos0 o , 0^ is the solar zenith angle, and 8 is 
the Earth-Sun distance correction factor. The narrow- 
band albedo is 

a 6 - o - C) x cs (9<, 9. ♦) P„ + Cx c (0„. 8. ♦) P s © 

where C is the total cloud amount, % cs and % c are the 
anisotropic reflectance correction factors for clear sky 
and clouds, respectively, from Minnis and Harrison 
(1984b). The mean dear-sky and cloud reflectances for 
the given scene are p Ci and p c , respectively. The respec- 
tive viewing zenith and relative azimuth angles are 0 
and<t>. 
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GOES-7 

During April 1994, GOES-7 was located in nearly 
the same position as GOES-6 had been during April 
1985. The 4-km VIS and IR data were combined into 
regional two-dimensional histograms for the 0.3° and 
0.5° grids described previously. The histograms were 
developed from data taken every half hour between 1330 
and 2330 UTC. The GOES-7 VIS counts D 1 were con- 
verted to reflectance by using the following equation: 

p 7 = (0.0126Z) 7 -4.0)/526.2 h o 8 (3) 

The calibration coefficients are based on an inter- 
calibration between GOES-7 VIS and the NOAA-11 
Advanced Very High Resolution Radiometer (AVHRR) 
channel 1. The procedure, described in appendix A, 
assumes that the AVHRR is calibrated. The equivalent 
blackbody temperatures were determined with the stan- 
dard GOES calibration. 

ERBS 

The Earth Radiation Budget Satellite (ERBS) mea- 
sured broadband regional shortwave (SW) and longwave 
(LW) radiances with a cross-track scanner having a nom- 
inal =31 x 47 km 2 field of view at nadir (Barkstrom and 
Smith 1986) as part of the Earth Radiation Budget Exper- 
iment (ERBE). The ERBS is in a 56° inclined orbit, 
which allows its equator-crossing time to precess all 
hours of the day in 36 days. The radiances were con- 
verted to SW and LW fluxes with the methods of Smith 
et al. (1986) and Wielicki and Green (1989) and the 
anisotropic models of Suttles et al. (1988, 1989) for 
clear, partly cloudy, mostly cloudy, and overcast pixels. 
Longwave fluxes within a given 2.5° latitude-longitude 
region, averaged to obtain the mean longwave flux Mj w 
for each ERBS overpass, are assumed to correspond to 
the local half hour at the center of the region. The SW 
albedo a w was computed from the SW flux and adjusted 
to the local half hour with the techniques of Brooks et al. 
(1986). The regional broadband quantities taken during 
April 1985 over the GOES-6 grid are used here. 

Atmospheric Profiles and Surface Data 

Vertical profiles of air temperature T a (p) and humid- 
ity from the National Meteorological Center gridded 
analyses at standard pressure p levels (surface; 925, 850, 
700, 500, 400, 300, 250, 200, 150, and 100 hPa; and 
tropopause) were used to analyze the April 1985 dataset. 
National Weather Service 12-hourly temperature and 
humidity soundings taken at standard levels during April 
1994 were interpolated to the mesoscale grid with a fast 
Barnes interpolation method (Hibbard and Wylie, 1985), 
as implemented by the University of Wisconsin 


Man-computer Interactive Data Analysis System 
(McIDAS). The interpolation technique uses all available 
soundings within the domain and the closest soundings 
outside the domain to eliminate edge effects. Hourly sur- 
face air temperatures T s from all reporting stations within 
the domain are interpolated to provide a default surface 
temperature for each grid box. 

Methodology 

The remote sensing of cloud properties involves the 
use of an idealized conceptual model of clouds to inter- 
pret the radiances emanating from real scenes. The con- 
ceptual model relies on the following assumptions. All 
clouds are plane-parallel entities that only occur within 
the boundaries of a satellite field of view or pixel. The 
clouds completely fill the pixels in which they occur and 
have a uniform distribution of particle sizes within the 
pixel. Each cloudy pixel contains a cloud at only one alti- 
tude and has a thickness prescribed by empirical formu- 
lae. The satellite sensor is correctly calibrated and the 
anisotropy of the clear-sky conditions is completely 
described with the assigned bidirectional reflectance 
model. The anisotropy of the cloud is completely 
described by the cloud bidirectional reflectance models 
used by ERBE for the SW data and by Minnis and 
Harrison (1984b) for the VIS data. The limb darkening of 
the IR radiances and LW fluxes follows the models 
developed by Minnis and Harrison (1984b) and the 
ERBE models, respectively. The details of the applied 
methodology are given in this section. Some short- 
comings due to these assumptions are mentioned in the 
“Discussion” section. 

Cloud Properties 

The methodology used here is designated the layer 
bispectral threshold method (LBTM). Cloud fraction Q, 
cloud-center temperature T ci , cloud-top temperature T ti , 
cloud optical depth T;, cloud-top altitude z ti , and cloud- 
center altitude z ci were computed for the following three 
height intervals: low (z c < 2 km), middle (2 < z c < 6 km), 
and high (z c >6 km). These computations were per- 
formed with the procedures described by Minnis, Heck, 
and Young (1993). The subscript i =1, 2,3 refers to the 
low, middle, and high layers, respectively. The total 
cloud amount C is the simple sum of the three layer 
amounts with no overlap, and the total cloud optical 
depth t, temperatures T t and T c , and altitudes z t and z c are 
weighted averages with C,- serving as the weights. The 
brightness temperature averages are performed with the 
equivalent radiance evaluated at 11.5 pm. The cloud- 
center temperature is defined as the equivalent radiating 
temperature of the cloud. For optically thin clouds, the 
cloud-center temperature generally corresponds to some 
temperature between the physical center and top of the 
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cloud. The cloud-center altitude z(T c ) is determined from 
a sounding. The cloud-center temperature approaches the 
cloud-top temperature as the cloud becomes optically 
thick. Thus, for optically thick clouds, T c = T t and z c = z t . 
A summary of the cloud analysis technique is given as 
follows. 

Clear-sky reflectance p CJ is determined at each hour 
for a given region with a minimum reflectance method 
(Minnis and Harrison, 1984a). These reflectances are 
assumed to be valid throughout the month; however, they 
may be adjusted for specific cases as in Minnis, Harrison, 
and Gibson (1987). During daylight hours, clear-sky 
temperature T cs corresponds to the mean IR radiance for 
all pixels having VIS counts D < D cs + AD and tempera- 
tures T>T[ im , where D cs corresponds to p M and 
AD = 10.4 + 1 .4 In (p 0 ). The term T[ im is the climatologi- 
cal minimum limit imposed on the clear-sky temperature 
(e.g., Minnis, Harrison, and Gibson 1987) to ensure that 
T cs is not unrealistically low. The clear-sky temperature 
is tested and altered as follows. The parameter is set 
equal to the smaller of the two values T s - 5 K or 
where 7/ iml = T S -AT SCS and A T scs is the mean hourly 
regional difference T cs - T s minus two standard devia- 
tions of that difference for a given local hour. Figure 2 
gives the mean values of T cs and T s for all scenes in the 
mesoscale domain having C < 0.05. The domain mean 
clear-sky temperature peaks near 1300 LT, while T s has a 
broad maximum approximately 2 hours later. Figure 3 
shows the differences between the mean values of T cs 
and T s minus two standard deviations of the difference 
(T cs - T s ) at a particular local time. The regression fit to 
the data shown in figure 3 is 

AT scs = - 19.5 + 3.2/1 — 0.144/t 2 (4) 

where h is the local time (LT) in hours and is referenced 
to 98°W. A second limit T / ;m2 = X s — 10 K is also used to 
constrain T cS . When a value of T cs is derived from the 
data, it is reset to equal T/ im if T cs < Ti iml . If no value of 
T cs can be determined (i.e., overcast conditions) and no 
temperatures are greater than 7/, ot2 or no VIS counts are 
lower than D cs + AD, then T s is substituted for T cs . If no 
value for T cs can be determined, but VIS counts are less 
than D cs + AD and observed temperatures are greater 
than T lim2> then T cs — T liml . 

The cloudy pixels are identified with a simple 
threshold approach that is based on the division of the 
VIS-1R histogram into clear and cloudy layers. (See 
fig. 4.) A pixel is cloudy when T < T cs - 5 K or 
D > D cs + AD. Cloudy pixels are grouped into low, mid- 
dle, and high categories with model calculations of tem- 
perature and reflectance variations for clouds at 2 and 
6 km. The calculations determine the p - T pairs that 
constitute the boundaries of each layer. These boundaries 


are denoted by the curves labeled P m , P h , and P p , which 
are shown in figure 4. The low-middle boundary P m is 
determined with a cloud composed of water droplets hav- 
ing an effective radius of 10 pm at the temperature T m . 
The middle-high boundary P h is computed with a theo- 
retical cirrostratus (CS) cloud at the temperature T h . The 
CS cloud comprises randomly oriented hexagonal ice 
columns having a range of sizes (Takano and Liou 1989). 
To constrain the interpretations, an upper boundary P p 
that corresponds to a CS cloud at the tropopause temper- 
ature T p is also computed. Pixels that are cold and darker 
than the P p boundary are defined as dark pixels. 

The model used to determine the boundaries and to 
interpret the reflectances for cloud optical depth is 

p = £p/(l-0 i = 1,5 (5) 

where £ is a regression normalization factor, and p,- are 
parameterizations of the multiple scattering and absorp- 
tion by the atmosphere, scattering by the cloud, and 
reflection by the surface (Minnis, Liou, and Takano 
1993). The reflectance parameterization is described 
briefly as follows. 

The visible-channel reflectance contributed by the 
cloud and the atmosphere above it is 

P'l “ f alPcl = f alPcl ^a) 

where p cl is the combined reflectance of the cloud and 
the Rayleigh scattering layer above it. The transmittance, 

t al = exp [-t^ (l/p 0 + 1/p.) ] 

where (X = cos 0, the Rayleigh optical depth above the 
cloud is and the ozone absorption optical depth for 
the VIS channel is x 0 . For this analysis, % Q is fixed at 
0.022, a value that corresponds to an ozone path length of 
0.32 cm-STP. The beam reflectance by the surface is 

p 2 = r c i t c T p, (6b) 

where the downward and upward cloud transmittances 
are 

t c l = exp [-( 1 -f D )t/p 0 ] 
and 

t c T = exp [-( 1 -f D ) t/(x] 

and f D is the fraction of the beam scattered in the forward 
direction because of diffraction and direct transmission 
through the droplet or crystal. The value of/ D is gener- 
ally greater than or equal to 0.5 at visible wavelengths. 
The fraction of radiation scattered from the forward 
direction, reflected by the surface, and transmitted 
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diffusely back through the cloud to space is approxi- 
mated as 

p 3 = aj\ - a cd )(\ - t c t - <x c ) (6c) 

where a c (x, |x 0 ) is the cloud albedo and a. cd is the diffuse 
cloud albedo. The fourth term 



accounts for the relative thickness of the Rayleigh layers 
above and below the cloud. The effects of the two 
Rayleigh layers are included by using the direct Rayleigh 
reflectance term for the bottom layer p R2 and the 
Rayleigh albedo for the top layer o. R \. The fifth term 

p 5 = fl 0 + a l( -hlVXs + a 2 a sd < 6e ) 

accounts for an overestimate in the surface contribution 
to the reflectance by p 2 for small cloud optical depths. 
The coefficients a,- depend on the microphysical model, 
a cs is the clear-sky albedo, and a sd is the diffuse clear- 
sky albedo. The denominator in equation (5) uses the fol- 
lowing parameter: 

C, = h 0 + hjln(T) +h 2 ct^ln(T) + bjCL sd (7) 

to minimize biases in the parameterization. The co- 
efficients bi also vary with the microphysical model. 

Cloud radiating temperature is assumed to be 

B(T) = (1 -e)B(T cs ) +eB(T c ) (8) 

where B is the Planck function evaluated at 11.5 pm. The 
effective emissivity is 

e = 1 - exp \_a (x/p) b ~\ (9) 

The coefficients a and b depend on the cloud micro- 
physics (see Minnis, Liou, and Takano 1993). 

The optical depth for the cloudy pixels is obtained by 
matching the observed reflectance to the reflectance 
parameterization through an iterative process. Optical 
depth cannot be computed for pixels as dark or darker 
than the clear scene, so they are averaged with brighter 
pixels until an optical depth can be computed. Cloud 
temperatures are computed using equations (8) and (9) 
with the observed IR temperature T cs and the derived 
value of X. If the value ofT c <T p -2 K and x/p < 5, then 
T c is reset to equal T p and the emissivity and optical 
depth are recomputed. Further details of these parameter- 
izations, the models used in them, and their application 
are found in Minnis, Heck, and Young (1993) and 
Minnis, Liou, and Takano (1993). 


Cloud-top temperature is computed with the cloud- 
top emissivity parameter. For optically thick clouds 
(x > 6), the cloud radiating center and physical top are 
assumed to be the same. Thus, the cloud-top emissivity 
E, = £. For optically thin clouds, £, depends on T c . If 
T c < 245 K, then 

e t = (- 0.00914T, + 2.966) £ (10) 

(See Minnis, Harrison, and Heck 1990). If T c > 280 K, 
then £ t - 0.99 e is assumed. For 245 K < T c < 280 K, lin- 
ear interpolation between the warm and cold cases, the 
following equation is used to determine the cloud-top 
emissivity: 

£ f = (0.00753 T c + 1.12) £ (11) 

For clouds having moderate optical depths (2 < x < 6), 
simple linear interpolation between the thick and thin 
estimates is used. "Hie cloud-top temperature is 

T t = IT 1 { [B(T) — (1 -£ t )B (TJ ] /£,} (12) 

where B~ l is the inverse Planck function. If T t < T p , then 
T t is reset to equal T p . Cloud-top altitude is taken from 
the sounding as z t = z(T t ). 

Cloud thickness was computed for clouds having 
T c < 245 K as 

A z = 7.2 - 0.024T C + 0.95 In (x) (13) 

(Smith et al. 1993) and for clouds with T c > 275 K as 

Az = 0.085 x 1/2 (14) 

(Minnis, Heck, et al. 1992). For cloud temperatures 
between 245 and 275 K, linear interpolation between 
equations (13) and (14) is used to compute Az. Several 
criteria are used to ensure that the cloud thickness is 
reasonable. If Az < 0.1 km, then Az is reset to equal 
0.1 km. If z t - Az < z s , then the cloud thickness is reset to 
Az-z t -z s + 0.1 km. If z t - Az> z c , then the cloud-top 
temperature is reset in the following manner: 
z t = z c + Az-0.1 km. These steps are repeated if 
necessary. 

Radiation Properties 

The SW albedos and top-of-the-atmosphere (TOA) 
LW fluxes for the ARM IOP are computed from the VIS 
albedos and IR equivalent blackbody temperatures, 
respectively, by using relationships based on regression 
analyses applied to the April 1985 GOES-6 and ERBS 
datasets. The ERBS data are matched to GOES-6 data 
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taken at the nearest UTC hour. The SW albedo is deter- 
mined as 

= a JVV (CLR) ( 1 - C) + a jw (CLD) C (15) 
where the clear-sky SW albedo is 


mesoscale region dataset (fig. 5(b)) is close to that deter- 
mined with all land regions in the GOES field of view. 
Combining the two fits as in equation (15) yields an 
overall uncertainty of 10.9 percent. 

The TOA LW flux is 


a ™( CLR ) = « 0 + a l a v( CLR ) +a 2 ( 16 > 

and the cloudy sky albedo is 

a iW ,(CLD) = b Q + b l a v (CLD) + & 2 a v (CLD ) 2 

+ h 3 ln(l/p 0 ) (17) 


M lw = c Q + c x M ir + c 2 M 2 ir + c 3 M ir In (RH) (18) 

where RH is the mean relative humidity for all layers 
above the altitude z(T a ) corresponding to T, as reported 
by Minnis, Harrison, and Young (1991). The narrowband 
flux is 


where (CLR) and Oy (CLD) are the clear and cloudy 
VIS albedos, respectively, and a t and b t are the regres- 
sion coefficients. The VIS albedos are restricted to val- 
ues between 0 and 1. If a v > 1, then its value is reset to 1 
before applying equation (17). The regression formula is 
unable to compute SW albedo greater than unity. These 
formulations yield approximately the same accuracy as 
the ray-matching techniques of Minnis and Harrison 
(1984b) because of the increased number of samples 
(Doelling et al. 1990). They also permit the determina- 
tion of regionally dependent coefficients. 

Preliminary studies of the regression analyses 
revealed that the clear-sky S W albedo can be more accu- 
rately determined if the regressions are performed by 
using matched clear-sky data and all data separately. 
Clear-sky data for the 2.5° regions are defined as those 
having C <0.15 from the HBTM and a clear-sky identifi- 
cation by the ERBE analysis. For the clear-sky case 
(eq. (16)), the quadratic term is not used because of the 
narrow range of data. Figure 5 shows the VIS-SW 
regression analyses for April 1985 with the clear-sky 
(fig. 5(a)) and all (fig. 5(b)) matched data with Q 0 < 78°. 
For figure 5(a), the multiple correlation coefficient 
R 2 = 0.87 and the standard error of the estimate (SE) is 
5.8 percent for data with 0 < 70°. The fit to the total 
dataset yields tf 2 = 0.91 and SE = 11.7 percent. The 
derived values for the clear-sky coefficients are 
a 0 = 0.0893, a\ = 0.5775, and a 2 = 0.0709. For the total 
dataset, only ERBS data taken for 0 < 45° were used to 
prevent the smearing effects of the larger ERBS fields of 
view at higher viewing angles. 1 The resulting coeffi- 
cients are %= 0.0588, b\- 0.8623, b 2 = -0.1190, and 
£>3 = 0.0624. Although the clear-sky curves are relatively 
flat, significant dependence on Q 0 is indicated by the 
separation of the plotted curves in both figures 5(a) 
and 5(b). The relationship derived with only the ARM 


1 Paper entitled “Reconstruction of Earth’s Radiation Field from 
Earth Radiation Budget Experiment (ERBE) Measurements” by 
Rajeeb Hazra, G. Louis Smith, and Stephen K. Park (pending 
publication). 


M ir = 2L ir 


c 2n r rc/2 

(0°) y. sin0 cos0 dQdty 

Jo J o 


= 6.18L iV (0°) 

(19) 

where the factor of 2 (in jxm) is used to account for the 
bandwidth (10.2 to 12.2 |tm) of the GOES sensor. The IR 
radiances 

II 

to 

3 

(20) 

are adjusted to the nadir view by 


M°°) = VM0) 

(21) 

where the limb-darkening function is 


y ir = 1, if 0 < 11° 



else 


y ir = 1.00067 + 0.03247 In (|Li) (22) 

This formula approximates the IR limb-darkening 
function used by Minnis and Harrison (1984a). Calibra- 
tion of the instrument accounts for the response function 
of the IR filter. 2 Thus, a flat sensor response is assumed 
when the equivalent blackbody temperatures are used. 
Equation (12) is an approximation of the narrowband 
flux developed to obtain the longwave flux from the cali- 
brated narrowband radiances. It was not intended to 
serve as a standard for comparison with calculations of 
window flux. The radiance corresponding to the equiva- 
lent blackbody temperature at the band center is probably 
a better standard for comparison. According to a personal 
communication with Menzel, the equivalent blackbody 
temperature was determined by using the band center to 
invert the measured radiance within ±0.1K of the mean 
temperature based on the radiance integrated over the 
entire band. 


Reference work on Prelaunch Study of VAS-D Performance 
that was performed by P. Menzel under contract NAS5-21965. 
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Figure 6 shows the matched IR-LW data and the 
resultant regression curves for various relative humidi- 
ties. These results are similar to those found for all land 
areas viewed by GOES during April 1985 (Minnis, 
Harrison, and Young 1991). The SE is 4 percent and R 2 
is 0.91 for the 1320 samples used in the correlation. Val- 
ues for the coefficients are c 0 = 64.39 W-m -2 , c\ = 6.57, 
c 2 - -0.0275 m 2 -W -1 , and C 3 = -0.322. These coeffi- 
cients were used in equation (18) to compute the LW 
fluxes from the April 1994 ARM IOP GOES-7 data. 

Results 

The methodology described in the preceding section 
was applied to the April 1994 GOES-7 data. The results 
have been archived for public distribution. Appendix B 
describes the format of the archived dataset and the 
means to access the data via the World Wide Web. 

Instantaneous Results 

To demonstrate the resolution of the cloud products, 
an example of the 1930 UTC, April 14, 1994 VIS imag- 
ery in figure 7 is overlaid with the derived values of T for 
a portion of the mesoscale grid. Its IR counterpart is 
shown in figure 8 overlaid with the cloud fraction in each 
region. These images contain fields of scattered thin cir- 
rus that are faint in some regions and thicker cirrus in the 
northwest and southeast comers of the image. Optical 
depths range from 0.1 to 8.7 in figure 7. The larger opti- 
cal depths are generally associated with the heaviest 
cloud cover (fig. 8 ). Figure 9 shows the entire mesoscale 
domain outlined in a VIS and IR image pair taken at 
1800 UTC, April 25, 1994. The primary feature includes 
several lines of thunderstorms passing northward through 
west central Oklahoma into Kansas. This squall line is 
preceded on the eastern side by a large mass of densely 
packed, low-level cumulus clouds and followed by thin 
wispy cirrus clouds in northwestern Texas. The thin cir- 
rus are evident in the IR image but barely discernible in 
the VIS photograph. 

The cloud properties derived from the images in fig- 
ure 9 are shown in figures 10 to 15. Total cloud amount 
(fig. 10 ) is nearly 100 percent over most of the domain. 
The clearer areas southwest and northeast of the thunder- 
storms are interspersed with thin midlevel and high 
clouds. The cirrus clouds in the Texas Panhandle were 
found primarily between 5 and 7 km. The lowest clouds 
are found in abundance over Arkansas and southwestern 
Missouri. These locations are consistent with the VIS 
imagery, which suggests that the cumulus clouds in those 
areas were dimmer and probably less developed verti- 
cally than the surrounding clouds. The surrounding 
clouds were placed in the midlevel category. Small frac- 
tions of low cloudiness are also observed along the 


periphery of the main high and midlevel cloud fields. 
Some small low-cloud amounts may result from pixels 
that are partially filled with midlevel or high clouds. 
Most midlevel cloud tops in the eastern region are 
located at ~3 km. 

Cloud optical depths (fig. 11) are generally very high 
over the domain, especially for the high clouds over 
Oklahoma and central Colorado. The mean cloud-center 
and cloud-top heights are shown in figure 12. The high- 
est clouds are found over the thunderstorms. The differ- 
ences between z t and z c are typically only 0.3 to 0.5 km. 
As expected, the thickest clouds (fig. 13) are found over 
Oklahoma and Colorado because of the cumulonimbus 
systems. Cloud thickness is estimated to be as great as 
6.5 km in some storms. All low clouds appear to be less 
than 1-km thick. Some of the densest midlevel clouds 
found between the high clouds are up to 3.5-km thick. 
The VIS and SW elear-sky and total albedos are shown 
in figure 14. The clear-sky VIS albedo varies from 0.10 
in central Arkansas to 0.18 in eastern New Mexico. 
Shortwave elear-sky albedos vary from less than 0.15 to 
over 0.20 in a pattern like that seen for the VIS albedo. 
The 1800-UTC clear-sky albedo near the central facility 
is =0.170. VIS total albedos are greater than the SW albe- 
dos except for the relatively clear regions. The total 
albedo variability is very similar to that for optical depth 
(fig. 1 1). Clear-sky IR and total temperatures are plotted 
with the corresponding LW fluxes in figure 15. The 
greatest Values of T cs and M !wcs are found in the clearest 
regions of the domain. Coldest cloud-tops with 
T c < 220 K are seen over the squalls in central 
Oklahoma. The LW fluxes vary from 120 to 305 W-m -2 . 

Average Cloud Properties 

The mean total and three-layer cloud amounts and 
optical depths averaged for a 2.5° region centered over 
the SCF are plotted in figures 16 and 17, respectively. In 
these figures, the results are shown for the 0.5° region 
including the SCF and the middle 0.3° box centered on 
the SCF. High cloudiness peaks during the late morning 
and late afternoon in the 2.5° box, while the midlevel 
clouds are relatively constant during the morning and 
diminish during the afternoon. Low clouds peak during 
the midaftemoon and gradually diminish toward the 
evening. Total cloud cover has a broad maximum cen- 
tered near 1330 LT. High clouds appear to have the 
greatest influence on the total cloud variability. In the 
0.3° box, the total cloudiness has a stronger early after- 
noon maximum. The high clouds show a broad maxi- 
mum from midmoming to late afternoon, while the 
midlevel clouds peak near 1300 LT. The low-cloud vari- 
ation is similar to that of the larger region but much nois- 
ier. Except for the noise in the 0.3° data, the cloud 
amounts for the two regions are similar. This similarity 
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•suggests that the surface measurements of clouds taken at 
the SCF may represent a GCM-scale gridbox for this area 
over the time scale of 1 month. 

Total optical depth (fig. 17) appears to be dominated 
by the high clouds at both scales. It peaks near 1330 LT 
then decreases rapidly an hour later, which suggests that 
the deepest daytime convection ends shortly after local 
noon. After the midday maximum in T, the high clouds 
comprise more thin cirrus blowoff than cumulonimbus 
tops. The midlevel clouds follow the mean total optical 
depths, and the low-level clouds are generally thinner 
than the higher clouds. An exception is the extremely 
high value at 1200 LT for the 0.3° region where the mean 
low-cloud optical depth reaches 34. This anomalous 
value resulted from a slight increase in the cloud-top 
temperature for a short time at midday during April 1 1 so 
that a few pixels were entered into the low-cloud cate- 
gory. A raining middle cloud layer with large optical 
depths covered the region at an altitude slightly higher 
than 2 km for most of the day. The hour-to-hour variabil- 
ity in the small region is absent in the 2.5° results. 

Figure 18 shows the mean half-hourly IOP total 
cloud amounts, heights, and optical depths averaged for 
the 2.5° and 0.3° regions in figure 17. The same quanti- 
ties are also shown for the 0.5° region that includes the 
SCF. The diurnal variations in the total parameters are 
similar in phase at all scales but differ in magnitude for 
the three regions. Hour-to-hour variability diminishes as 
the regional scale is decreased. The diurnal range in 
cloud height is =2 km in the smallest region; it drops to 
less than 1.5 km for the 0.5° region and less than 1.0 km 
for the 2.5° box. Similar decreases are evident in the 
mean total cloud optical depth and amounts. 

Figure 19 shows the distribution of observed cloud- 
top heights for the 0.5° region including the SCF. Clouds 
were observed most frequently between 2 and 5 km and 
between 9 and 13 km over the SCF. Few clouds were 
observed between 7 and 9 km. The corresponding cloud 
thickness distributions in figure 20 show that most of the 
low clouds were relatively thin with the mode thickness 
near 0.15 km. Midlevel clouds had a mode thickness of 
=0.20 km but had a much greater range than the low 
clouds. The high-cloud mode thickness is 0.5 km, but the 
range in thickness is about double that for the midlevel 
clouds. These results are similar to those for the 0.3° 
region centered on the SCF. 

The April IOP mean cloud properties for the SGP 
domain are shown in figures 21 to 24. Total cloud 
amount (fig. 21) varies from less than 40 percent in the 
southwestern comer to more than 70 percent in Louisi- 
ana, Colorado, and Nebraska. Low clouds appear to be a 
significant component only over western Oklahoma and 
the southeastern quadrant of the domain. Midlevel and 


high clouds dominate elsewhere, especially in the north- 
ern half of the grid. The cloud optical depth patterns 
(fig. 22) reveal the main centers of low-level conver- 
gence. Low, midlevel, and high cloud optical depths all 
reach maximum levels in southern Kansas and central 
Missouri. During this time period, convective storms 
such as those in figure 9 apparently form most often in 
the vicinity of the SCF and move northwestward through 
Kansas and Missouri. High clouds passing over the 
southwestern quadrant are primarily optically thin cirrus. 
Mean total cloud-center heights (fig. 23) range from 
=4 km along the southern and southeastern edges of the 
domain to almost 7 km over the Rocky Mountains. Other 
relative height maxima occur in a southwest-to-northeast 
line from the Texas Panhandle through the SCF and into 
Missouri. 

The mean low-cloud tops vary from 1.5 to more than 
2 km, while mean midlevel clouds are found between 3.5 
and 4.5 km. Although low clouds are only supposed to be 
found at altitudes below 2 km, some resetting of the 
cloud heights occurs after the initial analysis is per- 
formed. For example, the surface elevation (fig. 1) in 
some regions exceeds 2 km. Cloud height is determined 
with interpolated soundings. Thus, some interpolated 
values result in low-cloud altitudes over some of the 
higher elevations. When the cloud base (top-thickness) 
for a given cloud layer falls below the surface elevation, 
both cloud base and top are adjusted to levels at least 
0.1 km above the surface. Cloud height increases from 
southeast to northwest for low and midlevel clouds. 
Mean high-cloud heights, however, decrease from nearly 
11 km in the south to 8.5 km in the northern part of the 
grid. Figure 24 shows the variation in cloud-amount- 
weighted mean cloud thickness. Because of the heavy 
reliance of A z on optical depth, the patterns in figure 24 
are similar to those in figure 22. The thickest clouds 
occur in southern Kansas and central Missouri along the 
apparent mean path of the heavier thunderstorms. Mean 
mid-cloud thickness ranges from slightly less than 
0.6 km in the south to almost 2 km over the Colorado 
Rockies. Mean low-cloud thickness varies from 
=0.25 km to 0.35 km. 

Average Radiation Parameters 

Figure 25 shows the mean diurnal variation in the 
clear-sky, total, and effective cloud VIS and SW albedos 
for the 2.5° and 0.3° regions centered on the SCF. The 
mean SW effective cloud albedo is 

a = ra (LT) - a (CLR,LT) (1-C)1/C (23) 
swce L J 

Both VIS and SW clear-sky albedos increase with solar 
zenith angle (time from local noon). The slight diver- 
gence between the SW and VIS clear albedos during the 
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early morning and late afternoon arises from the 
dependence of equation (16) on the solar zenith angle. 
The total albedos increase only slightly with 0 O . Relative 
minima in the 0.3° data at 1030 and 1530 LT roughly 
correspond to the dips in the mean optical depth (fig. 18). 
Total VIS albedo is greater than the SW albedo around 
midday and less than the SW albedo at other times as a 
result of the 0 O dependence in equation (17). The effec- 
tive cloud albedo is generally higher during the morning 
than during the afternoon. For example, at 0900 LT a wce 
is 0.48 compared with 0.45 at 1500 LT for the same solar 
zenith angle. However, in figure 17, the mean optical 
depth at 1500 LT is greater than or equal to its 0900 LT 
counterpart. This apparent discrepancy results from the 
occurrence of two cumulonimbus events during the 
month at 1500 LT that produced x > 90. A comparison of 
the daily optical depths and effective cloud albedos at 
these two hours indicates that x < 3 for almost half the 
IOP days at 1500 LT compared with less than 25 percent 
of the days at 0900 LT. During most days, x (1500) < x 
(0900). TIius, the two thunderstorms increased the mean 
optical depth to values comparable to that at 0900 LT. 
However, the thunderstorms had a much smaller impact 
on the mean albedo because of the nonlinear relationship 
between albedo and optical depth. 

The mean clear LW fluxes and IR temperatures in 
figure 26 vary by 14 W-m -2 and 8 K, respectively, dur- 
ing the daylight period for both small and large regions 
centered on the SCF. Maximum clear-sky LW flux 
occurs between 1300 and 1330 LT. The maximum mean 
total IR temperatures and LW fluxes occur at 1530 LT 
when the mean middle and total cloud amounts diminish 
(fig. 16). The minimum at 0730 LT is driven by the 
clear-sky values, while the broad relative minimum 
between 1100 and 1430 LT in the 0.3° region results 
from the maxima in high, middle, and total Cloudiness. 
Much of this diurnal structure in the LW flux is 
smoothed out in the larger region. 

The mean domain clear-sky and total albedos in 
figure 27 reflect the patterns of vegetation and the com- 
binations of cloud amount, optical depth, and latitude, 
respectively. The clear-sky narrowband albedos for 
heavily forested areas are generally less than 0.13 with 
greater values in the north. Mixed forest, farmland, and 
tail-grass prairie areas have VIS clear-sky albedos 
between 0.13 and 0.15, and oq, (CLR) for mixed prairie 
and wheat-growing areas varies between 0.15 and 0.17. 
Clear-sky VIS albedo is greater than 0.17 for the high 
plains’ desert (z s > 1.2 km in fig. 1) and steppe regions. 
The clear-sky albedo also increases slightly with latitude 
because of increasing 0„. In figure 27, the clear-sky albe- 
dos are generally greater than their counterparts in fig- 
ure 14 because the latter were measured near local noon. 
Shortwave clear-sky albedos are greater than their VIS 


counterparts in figure 27 by 0.02 to 0.05. The values of 
a iW , which vary from 0.17 to 0.22, are consistent with 
the theoretical models presented by Briegleb et al, 
(1986). This range in values exceeds the ERBE range of 
0.15 to 0.21 (see fig. A2) by 0.01. Total albedo is at a 
maximum near 39°N and in central Nebraska because of 
(I) maxima in total cloud amount (fig. 21), (2) the rela- 
tively high optical depths (x > 9), and (3) the greater 
solar zenith angles. More cloudiness was observed in the 
southeastern comer, but x < 9 in those regions. The low- 
est total albedos are found in the southwest. The SW 
albedos are greater or less than the VIS albedos for oq, 
less or greater than =0.30 as expected from the VlS-to- 
SW conversion formula (fig. 5). 

Figure 28 shows the mean IR temperatures and LW 
fluxes for the domain. Clear-sky temperatures range 
from 281 K near the Rockies to 300 K in the southwest. 
The clear-sky LW flux follows a similar pattern varying 
from =257 to 285 W-m -2 . These values result from an 
average taken over all daylight hours and are much less 
than those shown in figure 15 for observations taken near 
local noon. The lowest values of mean total temperature 
and LW flux coincide with maxima in the high-cloud 
amounts (fig. 21). The clearest, warmest areas in the 
southwest have maximum fluxes of =270 W-m -2 . Inclu- 
sion of nighttime fluxes and temperatures would consid- 
erably reduce the values in figure 28. 

Discussion 

The assumptions used to interpret the data leave 
much room for uncertainty in the derived products. 
Quantification of all uncertainties is not possible. How- 
ever, the sources of uncertainty are discussed here and 
estimates of their impact on a particular cloud property 
are noted when possible. 

Cloud Amount 

Clouds typically do not completely fill all pixels, and 
clear pixels are not always entirely free of clouds. The 
VIS and IR thresholds are set to some value above the 
expected clear-sky values. Thus, some pixels with small 
amounts of cloudiness are classified as clear and some 
pixels with clear portions are designated as cloudy. 
These pixels hopefiilly offset each other to yield the 
proper cloud amount for the entire area. Whenever only a 
few scattered clouds, very thin clouds, or a few clear 
holes are in a scene, probably no compensatory effects 
will occur due to the threshold selection. Wielicki and 
Parker (1992) performed a study of effects of pixel reso- 
lution on cloud amount using several different analysis 
algorithms. Although different thresholds and pixel aver- 
aging techniques are used, the LBTM is most similar to 
the ISCCP method highlighted in the Wielicki and Parker 


10 



.(1992) study. For the 24 Landsat scenes examined by 
Wielicki and Parker (1992), the ISCCP method using 
4-km resolution data overestimated cloud cover for vari- 
ous broken water-cloud scenes by 0.025 and under- 
estimated cirrus by 0.08. Similar values may be expected 
for the LBTM in like conditions. The Landsat study was 
specifically focused on broken and scattered cloud fields 
and very thin cirrus clouds. It was also restricted to the 
Landsat’ s near-nadir viewing angle. Over the course of a 
month, both thin and thick stratiform and cumuliform 
clouds will pass over a given area reducing the influence 
of broken cloud scene errors. The GOES viewing angle 
to the SGP locale is =44°. Thus, the overall effect of 
uncertainties due to resolution and pixel-filling are still 
undetermined. Evaluation of higher resolution data and 
comparisons with other observations are needed to assess 
the errors in cloud amount. 

Comparisons were performed with visual observa- 
tions from Whiteman Air Force Base (WAFB) at 
38.73°N, 93.55°W and Tinker Air Force Base (TAFB) at 
35.42°N, 97.38°W to examine the uncertainties in total 
cloud cover. These datasets were used because they 
included total sky cover estimates and comments about 
sky conditions. The hourly surface observations (taken in 
tenths of sky cover converted to percentages) C s f c were 
compared with values of C from the 0.5° region centered 
at 38.75°N, 93.75°W, from a 1.0° region centered on 
WAFB, and from the 1.0° region (average of two 0.5° 
regions: 0.5° longitude by 0.5° latitude) immediately sur- 
rounding TAFB. From the surface, the viewing radius for 
high clouds is =40 km and much smaller for low clouds. 
Figure 29 shows the mean LBTM cloud amounts found 
for each tenth of sky cover. The number of observations 
in each tenth is shown at the top of each plot. The LBTM 
data were linearly interpolated to fill in a few missing 
hours. Nearly complete agreement was obtained for all 
clear and overcast skies, as shown from the surface. Over 
WAFB (fig. 29(a)), two cases occurred for which the sur- 
face reported clear and the LBTM retrieved a few percent 
low clouds. These clouds are evident in the northwest 
comer of the box in the imagery and could not have been 
viewed by an observer at the center of the box. Although 
many cases occurred for which the surface reported less 
than 100-percent cloudiness and the satellite based 
cloudiness was 100 percent, no cases occurred for the 
opposite. Thus, the LBTM as formulated here generally 
will not produce cloud observations when there are no 
clouds and will not produce clear spots when none exists. 
Exceptions may arise over snow surfaces or in shadows. 

Figure 29(a) clearly shows that the LBTM tends to 
overestimate cloud fraction relative to sky cover for 
larger cloud amounts and underestimate for smaller frac- 
tions of sky cover. Surface observers record a minimum 


of 10-percent sky cover if any clouds are observed and a 
maximum of 90 percent if the clouds have any holes. 
Thus, some relative under or overestimation by the satel- 
lite is expected for the 10- and 90-percent sky cover cate- 
gories. However, the underestimation of sky cover for 
the range where 10 percent < Csfc < 60 percent is prima- 
rily the result of missing thin cirrus or some scattered 
small, low cumulus clouds. A detailed examination of the 
data over WAFB revealed that on some occasions these 
clouds did not produce a signal large enough to be 
detected in either the VIS or IR channels. In several 
cases, even extreme contrast enhancement failed to 
reveal the clouds in the satellite data. Detection of such 
clouds with an automated scheme involves a greater risk 
of false cloud identifications. The overall impact of miss- 
ing these extremely tenuous clouds on the total cloud 
fraction or on other aspects of the water or radiative bal- 
ance is minimal because of their relatively small optical 
depth, areal coverage, and water content. More cloud 
cover was found in the 1.0° region than in the 0.5° 
region. Sometimes the clouds viewed by the surface 
observer occurred to the east of WAFB so that they were 
missed in the 0.5° region. This difference is evident in 
the relative number of partly cloudy (10 to 40 percent) 
surface reports that were determined as clear by the 
LBTM. For the partly cloudy surface reports, no clouds 
were found in 41 and 16 percent of the satellite retrievals 
for the 0.5° and 1.0° regions, respectively. Thus, the 
comparison with the 1.0° region is probably a better 
guide to the relative accuracy of the LBTM. 

Much better agreement exists between the LBTM 
and the TAFB data (fig. 29(b)). Although the trends are 
similar to the WAFB results, the differences for 
Csfc < 60 percent are much smaller. Fewer overcast and 
clear cases and 15 percent more broken and scattered 
cloud conditions occurred at TAFB than at WAFB. 
Whether the differences in the two satellite-surface com- 
parisons are due to the particular observers at the two 
locations, the clouds, or the angular conditions is un- 
certain. Considerable evidence and logic exists for 
expecting some systematic differences between surface- 
observed and satellite-derived cloudiness. Henderson- 
Sellers and McGuffie (1990) explored this idea with all- 
sky camera photographs to simulate sky cover and Earth 
cover (satellite viewed cloudiness). They derived a rela- 
tionship in oktas that converts to 

C = 0.0459 + 0.0694C^ /c 

+ 0.0219C* /c - 0.0001 25 C^ /c (24) 

when C is expressed in percent. Equation (24) was used 
to compute the expected satellite-observed cloudiness 
and compared to the averaged decile cloud amounts from 
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WAFB and TAFB (fig. 30). Except for the disagree- 
ments at 40 and 80 percent, the empirical model is 
consistent with the LBTM-surface comparison. Thus, 
viewing perspective may explain much of the difference 
between the surface and satellite observations. 

Figure 31 presents the mean hourly cloud amounts 
for all times when there are observations for both the sat- 
ellite and surface sites. Except for the first three morning 
hours, the mean LBTM cloudiness for the 1.0° region is 
within ±2 percent of the WAFB averages (fig. 31(a)), 
The differences in the mean morning cloud amounts are 
due to the thin cirrus problem noted earlier. The agree- 
ment with the 0.5° region is not quite as good with 
8-percent differences at 1600 and 2300 UTC. Some dif- 
ferences between the surface and 0.5° region are due to 
the mismatched fields of view as noted earlier. All three 
curves are well matched during most of the afternoon. 
Overall, the satellite estimates for the 0.5° and 1.0° 
regions are 2.9 and 1.9 percent less than the WAFB sky 
cover, respectively. The hourly means have rms differ- 
ences of 4.3 and 2.9 percent, respectively. The average 
difference for the 1.0° region is less than half the mean 
difference found between International Satellite Cloud 
Climatology Project (ISCCP) daytime land cloudiness 
and surface observations (Rossow, Walker, and Garder 
1993). Slightly better correspondence exists between the 
TAFB observations and the LBTM results for the sur- 
rounding 1.0° x 0.5° region in figure 31(b). Here, the 
LBTM mean cloudiness is only 0.5 percent less than the 
surface. The hourly mean rms difference is 2.3 percent. 

Figure 32 shows the mean hourly, surface-observed 
cloudiness for the SCF compared with the averaged 
0.3° domain LBTM results (fig. 32(a)) and for Vance 
Air Force Base (VAFB) at 36.33°N, 97.92°W compared 
with two 0.5° grid boxes straddling VAFB (fig. 32(b)). 
VAFB is only =31 km from the SCF. Although the 
LBTM mean cloudiness exceeds the SCF data by 
3.1 percent, it is 0.9 percent less than the VAFB results. 
The hourly mean rms differences are 4.6 and 4.1 percent, 
respectively. The VAFB results are consistent with the 
other Air Force comparisons. The overestimate of the 
satellite cloud amount compared with the SCF may be 
due to differences in the surface observer’s reporting. 
The SCF observers report fractional cloudiness for each 
sky quadrant. These quantities were then averaged to 
obtain the total sky coverage. Also, the Air Force observ- 
ers may have been trained differently than the SCF 
observers. The SCF mean cloud amount is 8.6 percent 
less than the VAFB, which is almost twice the difference 
between the corresponding LBTM results. Overall, the 
mean LBTM cloudiness is 1.1 percent less than the mean 
for the three Air Force bases. Inclusion of the SCF results 
yields an average difference of -0.1 between the LBTM 


and surface results. The corresponding hourly mean rms 
differences are 3.2 and 3.6 percent, respectively. 

Given the assumption that the surface observer can- 
not determine cloud amount to better than ±15 percent 
for a given observation in scattered or broken cloud con- 
ditions, the uncertainty in the hourly mean for the 27-day 
period would be =3 percent. Thus, the LBTM hourly 
means would be within the uncertainty of the surface- 
observed hourly averages. However, the uncertainty in 
the surface-observed cloud fraction is probably less than 
15 percent because of the frequent occurrence of over- 
cast and clear skies, which are generally accurately 
observed. Perfect agreement with surface observations is 
also not expected because of both the subjective nature of 
the surface data and the cloudbase-height-dependent sur- 
face area that corresponds to the sky cover reports. 

Although the comparisons with the surface observa- 
tions are encouraging, they represent only four regions in 
the domain. Other factors such as snow cover, which 
occurred over some regions during the month, are mis- 
interpreted as clouds. If the snow is thick, it may be mis- 
taken for an optically dense cloud. The brighter surface 
in the western part of the domain may exacerbate the 
identification of thin cirrus clouds because of diminished 
contrast over barren surfaces. Until more effective mea- 
sures are developed for identifying snow and thin clouds, 
users of these data products must take such limitations 
into account. 

Cloud Height and Layering 

Cloud-center height errors for the LBTM were 
examined by Minnis, Heck, and Young (1993) for cirrus 
clouds over Wisconsin. The cirrus heights were within 
±1.3 km of the values derived from lidar images. With 
uncertainties in the lidar-derived altitudes taken into 
account, the uncertainty in the cirrus altitudes for clouds 
having t < 5 is about ±1.0 km. Similar results were found 
by Minnis, Young, et al. (1992) and Smith et al. (1993). 
All clouds having T c > 253 K are assumed to consist of 
water droplets. If the cloud is thin and consists of ice 
crystals, the altitude is underestimated by a factor that 
depends on the angles and optical depth. For thicker 
clouds, the errors in z c and z t are much smaller because 
the observed temperature of the cloud is close to the 
actual temperature. Thus, the vertical placement of the 
cloud will be somewhat dependent on the accuracy of the 
temperature profile. For single-layer clouds, the typical 
retrieved cloud-top or center height is estimated to be 
accurate to less than 1 km and 0.5 km for thin and thick 
clouds, respectively. For thick clouds over thick clouds, 
the errors are the same for the single-layer case. If a thin 
cloud overlays a thick cloud, the cloud-top altitude is 
underestimated because the algorithm interprets the VIS 
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data as an optically thick cloud. Thus, the derived alti- 
tude falls between the lower cloud and upper cloud. The 
error depends on the height difference between the two 
clouds, the upper cloud optical depth, and the viewing 
angle. Some validation of the cloud heights may be pos- 
sible by using the surface data taken during the April 
IOP. 

The WAFB observations are useful for examining 
how the layering of the clouds differs depending on the 
viewpoint. Three levels of cloudiness were recorded as 
scattered, broken, or overcast for the first three layers 
observable from the surface. Estimates of cloudbase 
height were used to assign the clouds to a layer compara- 
ble to the LBTM estimates. Cloud fraction at each level 
was assigned a value of 0, 33, 67, and 100 percent 
depending on its cloudiness categories. Using only the 
synoptic hours 1500, 1800, and 2100 UTC, multilayer 
and single-layer clouds were observed 42 and 36 percent 
of the time, respectively. Although the satellite detected 
clouds at more than one level 58 percent of the time, 
single-layer clouds comprised only 23 percent of the sat- 
ellite results. Low, middle, and high cloud amouhts from 
the surface were 31.2, 23.1, and 35.1 percent, respec- 
tively. The corresponding LBTM estimates for the 1.0° 
region (e.g., fig. 30) are 5.9, 24.8, and 31.3 percent. The 
low-cloud differences are primarily due to upperlevel 
clouds obscuring the low clouds from the satellite per- 
spective. The greater frequency of multilayer clouds in 
the satellite results arise from several factors. Variable 
high-cloud emissivity over a lower level cloud will yield 
a variable cloud height because the data are interpreted as 
being optically thick. The layering of the clouds above 
low overcast cannot be determined from the surface but 
can be estimated from the satellite. When a cloud’s alti- 
tude overlaps two different levels, two cloud layers are 
identified. Some partially high-cloud-filled pixels may 
be interpreted as lower clouds. These various factors 
make precise determination of cloud layering from either 
perspective difficult. 

To examine some differences in viewing perspec- 
tive, all cases containing low clouds were removed when 
one of the upper layers was reported as overcast or when 
rain was noted (assumes the clouds are thick enough to 
extend into the higher layers). The low-cloud fraction 
from the remaining surface observations reduces to 
8.1 percent, a value close to the satellite result. Middle 
clouds are eliminated when the highest level is reported 
as overcast. The middle-cloud amount is only given half 
weight when the low or middle levels are overcast and a 
thunderstorm is reported. If low overcast occurs with 
rain, the middle clouds are assumed to be overcast. This 
correction yields an estimate of 24-percent midlevel 
clouds from the satellite. Finally, the observed high- 
cloud fraction is altered by assuming that 50 percent of 


the sky is covered by high clouds when a thunderstorm is 
reported and lower clouds obscure the high clouds. The 
result is a net increase in high-cloud amount to 38 per- 
cent. The sum of these “adjusted” cloud amounts is 
70 percent, a value slightly greater than the mean sky 
cover. Such a correction is subject to a number of errors 
because of the assumptions made in the formulation. 
Nonetheless, the results are consistent with the satellite- 
derived cloud amounts in each layer. As previously 
noted, the greater amount of high cloudiness estimated 
from the surface data is probably due mostly to the thin 
cirrus clouds missed by the satellite analysis. While this 
comparison provides a somewhat quantitative verifica- 
tion of the layer cloud analyses, it clearly demonstrates 
the limitations in understanding the cloud field when 
only one view is available. 

Cloud Optical Depth and Thickness 

Cloud optical depth uncertainties are more difficult 
to evaluate. Some errors are expected to occur from the 
use of a parameterization instead of a fully detailed radia- 
tive transfer model. The parameterization errors are 
given by Minnis, Liou, and Takano (1993) and are less 
than 10 percent for the types of backgrounds in the SGP 
domain and for the angles observed in this dataset. The 
three-dimensional characteristics of real clouds also 
affect the optical depths retrieved with the plane-parallel 
model. The magnitude of this error is unknown. Minnis, 
Heck, and Young (1993) compared values of T derived 
by using simultaneous GOES and AVHRR data taken 
over inhomogeneous cirrus clouds from different angles. 
They found that the mean optical depths derived from the 
AVHRR angles were 9 percent greater than those 
derived from the GOES angles with the CS model. The 
rms difference is 45 percent. That comparison provides 
an overall estimate of the uncertainty in T for single-layer 
cirrus clouds. A similar study for water droplet clouds 
has not been conducted yet. The errors for broken water 
clouds are probably close to those found for the cirrus 
clouds because the cirrus clouds in the Minnis, Heck, and 
Young (1993) study were generally very inhomoge- 
neous. The errors are likely to be smaller over thicker 
stratiform clouds. If an ice cloud has T c > 253 K, its opti- 
cal depth is usually overestimated because of the water- 
droplet model used in the analysis for midlevel clouds. 
An examination of the ISCCP data for April 1984-1991 
indicated that cloud optical depth in the SGP locale is 
highly variable interannually. The optical depths in fig- 
ure 18 for the 2.5 region are within one standard devia- 
tion of the mean April 1984—1991 values for the closest 
ISCCP region. 

Cloud thickness is a remote-sensing product that has 
received minimal attention. The formulae used to com- 
pute cloud thickness are limited in scope and have not 
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yet been tested with independent datasets. The rms 
regression errors for equations (12) and (13) are ±0.92 
and 0.06 km, respectively, and were found using only 
single-layer cloud systems. It is not known if these errors 
apply to the data analyzed here. The impact of cloud 
overlap is also uncertain. Correlative measurements from 
instruments such as radar are needed to evaluate the 
derived cloud depths. 

The distributions of middle- and high-cloud thick- 
nesses in figure 20 are somewhat similar to those 
reported by Carswell et al. (1995) for lidar observations 
taken at Toronto, Canada. Although the Carswell et al. 
(1995) data were taken in a different area and period and 
only represent relatively thin clouds, the data provide an 
example of true thickness variability. The Toronto 
midlevel clouds have a thickness frequency that peaks 
between 0.27 and 0.50 km and decreases with increasing 
thickness to values greater than 3.0 km. The Toronto 
high clouds have a broad maximum in thickness between 
0.20 and 1.50 km with decreasing occurrences up to 
3.0 km. Those distributions are similar to the histograms 
in figure 20 without the thickest clouds. The tops of 
thicker clouds, however, could not be detected with the 
lidar. Maximum low-cloud thicknesses over Toronto 
occurred between 0.4 and 0.6 km. Those maximum val- 
ues are more than twice the values found in figure 20 for 
low clouds. Some difference may be caused by the view- 
ing perspective. Surface-observers determine low cloudi- 
ness from the base altitude, while the satellite uses the 
top altitude. Thus, some low clouds over Toronto could 
be designated as middle clouds by the satellite. More 
likely, the formula for estimating low-cloud thickness 
may underestimate cloud thickness over land because it 
is based on measurements of marine stratocumulus 
clouds. The comparison with the Toronto results does not 
verify the thicknesses derived in this study. It indicates, 
however, that satellite-derived values for the midlevel 
and high clouds are realistic and that a different approach 
may be needed for the lower clouds. Radar data taken 
during ARM IOP’s and other experiments will be useful 
for improving cloud thickness estimates. 

Radiation Properties 

Some clouds can act as plane-parallel sheets, espe- 
cially stratiform and optically thick clouds with small 
vertical aspect ratios (thickness/width). Cumuliform 
clouds and many cirrus clouds have complex shapes and 
may have relatively large (>0.5) aspect ratios. Thus, the 
radiances reflected by these latter cloud types may differ 
substantially from those by plane-parallel clouds. This 
effect is readily seen for cumulonimbus clouds and the 
edges of thick stratus clouds when the Sun is off the 
zenith and illuminates the sides of the clouds. The cloud 
appears brighter than normal when viewed from the solar 


side and darker than normal when viewed from the 
shaded side. An extreme example was observed at 2130 
UTC, April 10 in the 0.3° region centered over the cen- 
tral facility. The edge of a large cloud was illuminated, 
which resulted in a VIS albedo of 1.03. This value is 
unrealistic, especially since some visible radiation is 
absorbed by ozone. Such effects can be minimized by 
averaging over larger scales that include the tops as well 
as the sides of the clouds. Because the geostationary sat- 
ellite always views a given area at a constant angle, 
accounting for shading is not always possible, especially 
in extratropical latitudes where the relative azimuth 
angles are usually greater than 90°. The anisotropic 
reflectance models are based on observations and implic- 
itly account for some nonplane-parallel effects. How- 
ever, the models cannot account for the extreme 
variations from the mean that are observed in some 
cases. 

The anisotropic corrections for clear land are also 
general empirical models that may not be totally applica- 
ble to the particular surfaces and clear atmospheres in 
this dataset. However, the results indicate that the correc- 
tions are relatively accurate. Figure A3 (see appendix A) 
shows that the clear-sky reflectance is lower in the morn- 
ing than during the afternoon. Except for cases of frost, 
dew, or haze that are diumally dependent, the clear-sky 
albedo is expected to vary only with the solar zenith 
angle. Thus, the albedo should be symmetrical about 
local noon. Also, the surface and clear-sky albedos typi- 
cally increase with solar zenith angle, which results in 
higher values in the morning and evening and a mini- 
mum value near noon. The VIS clear-sky albedos in fig- 
ure 25 are almost symmetrical about the minimum at 
local noon. (The symmetry is evaluated in the next para- 
graph.) These albedos are based on clear-sky reflectances 
that are very similar to those in figure A3. Thus, the cor- 
rection factors have at least resulted in albedos that are 
consistent with the typical behavior. Simultaneous views 
from other satellites or from aircraft are needed to quanti- 
tatively assess anisotropic correction uncertainties for 
both cloudy and clear conditions. 

Errors in the LW fluxes and SW albedos due to 
the conversion process from narrowband-to-broadband 
fluxes are represented by the standard errors of the esti- 
mate given previously. Further comparisons with ERBE 
data, however, are needed to evaluate the data after the 
conversion formulae are applied. Clear-sky albedos 
derived from April 1985 ERBS data over the regions 
straddling the SCF are plotted as a function of \i a in 
figure 33 with the mean 1994 SW clear-sky albedos for 
the same area derived from the GOES-7 data. ERBS data 
were used only when more than 50 percent of the ERBS 
pixels were identified as clear. This comparison shows 
that the GOES-derived albedos agree well with the 


14 



.earlier ERBE observations at low solar zenith angles but 
diverge from the ERBE values at higher angles (lower 
*t 0 ). Both the GOES and the ERBE results can be fitted 
with a quadratic curve. The respective regression fits for 
the observed range of \x 0 are 

a sw (ERBE) = 0.341 - O.358p o + 0.201*1* 
and 

a sw (GOES) = 0.363 -0.356*t o + 0.161*1* 

The rms differences in the mean clear-sky albedos inte- 
grated over the observed solar zenith angle range is 
5.6 percent. Given that the 1-c interannual variability in 
clear-sky albedo is 5 percent over the domain (see appen- 
dix A), some difference may be due to interannual vari- 
ability. Some difference is probably due to errors in the 
anisotropic correction. Note also that the sampling of the 
ERBE data is relatively sparse possibly as a result of die 
ERBE scene selection scheme (Wielicki and Green 
1989), which rejects albedos greater than the preordained 
values. 

Uncertainties in the VIS channel calibration further 
influence the derived broadband albedos. Although the 
VIS calibration was shown to be reasonable, the accu- 
racy can be improved by using other methods to calibrate 
the GOES VIS channel. A bias in the albedos is likely to 
result from a calibration error. Indirect techniques such 
as those used by Whitlock, LeCroy, and Wheeler (1994) 
or direct methods that align well-calibrated, airborne, 
small field-of-view radiometers with the satellite view 
yield more reliable calibrations and a more accurate 
assessment of the uncertainties in the observed radiances 
than the approach used here. The IR-LW conversion is 
applied with the assumption that the GOES IR calibra- 
tions in April 1985 and 1994 are identical. Any differ- 
ence in the thermal channel calibrations will cause a bias 
in the derived LW fluxes. 


Concluding Remarks 

Version 1.0.0 of the Atmospheric Radiation Mea- 
surement (ARM) satellite data analysis program has been 
described and applied to Geostationary Operational 
Environmental Satellite (GOES) data taken during April 
1994. This initial step in the development of a compre- 
hensive analysis program yields reasonably accurate 
cloud properties given the limitations of field of view and 
spectral coverage of the satellite instrument. A technique 
for calibrating the GOES VIS (visible) channel has also 
been described and validated with dear-scene re- 
flectances. This cloud retrieval program is the first to 
take advantage of the extensive array of surface and 
sounding observations because it has been developed on 
the Man-computer Interactive Data Analysis System 
(McIDAS). The McIDAS has an extensive, near-real- 
time set of meteorological data. The availability of co- 
incident surface observations over much of the domain 
permits the application of more accurate constraints on 
the dear-cloud thresholds than previously possible. The 
results of this analysis provide the highest resolution 
satellite-derived cloud product available for ARM stud- 
ies and fill a gap between surface measurements and gen- 
eral circulation models (GCM). These data can be 
accessed on the Internet via the World Wide Web at the 
following uniform resource locator: 

http: //albedo . larc .nasa . gov: 1123/arm.html 

Future versions of this analysis technique will diminish 
some shortcomings of the present version by incorporat- 
ing additional spectral channels and provide a greater 
array of cloud products. New methods for incorporating 
more surface observations will also be explored. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
July 7, 1995 
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Appendix A 


Calibration of GOES-7 Visible Channel Using 
NOAA-11 AVHRR 

The GOES VISSR (Geostationary Operational Envi- 
ronmental Satellite Visible Infrared Spin Scan Radiome- 
ter) visible channel comprises eight sensors that have a 
nominal resolution of 1 km at nadir. The 4-km and 8-km 
resolutions used in the cloud and radiation analyses are 
obtained by computing the root mean square of die digi- 
tal counts for sets of 4x4 or 8 x 8 line-element arrays of 
contiguous data. The gain of each sensor may be periodi- 
cally adjusted to account for degradation, striping, and 
seasonal variations in the insolation at the top of the 
atmosphere. Thus, the overall calibration gain of the visi- 
ble data is subject to considerable variability. The infra- 
red (IR) sensor is nominally an elongated field of view 
with a nominal resolution of 4 km in the scanning (east- 
west) direction and 8 km in the line (north-south) direc- 
tion. The 4-km and 8-km resolutions are obtained by 
duplicating scan lines or averaging scan elements, 
respectively. An onboard thermal blackbody, which is 
used to monitor the IR channel gain, keeps the IR sensor 
calibrated. Except for the eclipse periods near the equi- 
noxes, the VISSR IR channel yields relatively stable val- 
ues of equivalent blackbody temperature. Therefore, the 
visible channel has the most acute need for calibration. 

The most recent calibration of the GOES-7 visible 
channel was performed during December 1991 by 
Whitlock, LeCroy, and Wheeler (1994) using a surface- 
based indirect method. Given the variability of the 
GOES-7 gain, simply using the 1991 value is not feasi- 
ble. To achieve a more applicable gain for April 1994, 
the NOAA-11 Advanced Very High Resolution Radio- 
meter (AVHRR) is used to perform an intercalibration 
with GOES-7. The NOAA-11 is in a Sun synchronous 
orbit with an equatorial crossing time of ~0400 and 1600 
LT during April 1994. Its visible sensor (channel 1) is 
not adjusted; however, its response may degrade with 
time. The NOAA-11 AVHRR channel 1 was calibrated 
several times by various methods between September 
1988 and December 1991. This time series of calibra- 
tions (Whitlock, LeCroy, and Wheeler 1994) was used to 
compute a correction for the sensor degradation that is 
assumed to have continued at the same rate until its fail- 
ure during September 1994. The AVHRR gain is 

g A = 0.567 4 - 0.0000352^ (Al) . 


where d s i is the number of days since launch on Septem- 
ber 24, 1988. The reflectance in percent for the Sun at 
lAUis 


p A = 5^ ( ^A- 4Q -°) CA2) 

where D A is the 10-bit AVHRR visible count. The 
GOES-equivalent radiance is 

L a = 5.269p A 8 (A3) 


The value of L A is then adjusted to correspond to a GOES 
view by correcting for the differences in the viewing and 
illumination angles as follows: 


v= 


x(e o ,e, V ) G n oG 

x(e o ,e, ¥ ) A n oA - 


(A4) 


where % is the mean anisotropic reflectance correction 
factor and the subscripts A and G refer to the GOES and 
AVHRR, respectively. The values of % are determined by 
weighting the clear-sky and cloudy factors by the clear 
and cloudy portions of the scene as determined from the 
GOES data. 

The April 1994 GOES-7 and NOAA-11 AVHRR 
data were matched on the 0.5° grid for the Atmospheric 
Radiation Measurement (ARM) mesoscale regime. Val- 
ues for L g were computed with equation (A4) and the 
average AVHRR channel-1 10-bit counts. These data 
were regressed against the root mean square GOES 
counts for the same regions taken within 15 min of the 
AVHRR overpass. The initial, unconstrained regression 
that was obtained by using all the data produced the dot- 
ted line in figure Al. The initial fit has R 2 - 0.93 and 
intercepts the X-axis at 417. During April 1994, the 
GOES visible offset count (4.7) was determined by view- 
ing space with the satellite. Thus, the regression fit 
should cross the X-axis at 22. The data were regressed 
again forcing the offset to the observed value. The results 
are shown in figure Al with the dashed line. Relatively 
poor sampling and considerable scatter occurred in the 
data for D 1 > 10000. The brighter points, however, have 
a significant effect on the regression line. The anisotropic 
correction for the clear land is expected to be more reli- 
able than the correction for clouds (brighter scenes), 
which will result in less scatter for the darker scenes. To 
minimize the effects of the unbalanced sampling and 
cloud anisotropy, only those data having D G <100 were 
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used in the regression. The resulting formula; shown as 
the solid line in figure Al, is 

L g - 0.0127Dg - 0.2804 (A5) 

This relationship, based on 300 matched data points, is 
used for the April 1994 ARM intensive observing period 
(IOP). 

To ensure that this approach is reasonably accurate, 
the resulting narrowband VIS clear-sky reflectances were 
compared with those determined over the same area dur- 
ing April 1985. Over any given land area, except in vari- 
able snow conditions, the clear-sky albedo should remain 
relatively constant from 1 year to the next. This interan- 
nual consistency is evident in the April clear-sky SW 
albedos computed from the ERBE data over the ARM 
Southern Great Plains (SGP) domain for 1985 to 1989. 
Figure A2 shows a scatterplot of the annual clear-sky 
SW albedos with the April 1985 to 1989 means for the 
2.5° regions within the SGP domain. The smaller albedos 
(v. sw <0.19) primarily correspond to the regions east of 
the central facility, and the others are west of the site. 
The variability in a sw is greater for the western regions 
presumably because of interannual differences in snow 
cover, soil moisture, and vegetation growth. The surface 
albedos in the eastern areas are probably less sensitive to 
soil moisture because of greater forest coverage. The 
overall interannual variability for the domain is ±0.011 
or =5 percent as measured at the two standard deviation 
a levels. Because the VIS albedo is directly proportional 
to the shortwave (SW) albedo, the interannual variability 
in <x w should be approximately the same for the visible 
(VIS) data. Thus, the clear-sky reflectance or albedo 
should not differ from 1 year to the next by more than 
10 percent (2-a level is 5 percent). 

Two datasets, the April 1985 GOES-6 data analyzed 
with the hybrid bispectral threshold method (HBTM) and 
by the International Satellite Cloud Climatology Project 
(ISCCP) by Rossow and Schiffer (1991) are used to eval- 
uate the April 1994 results. The ISCCP methodology 


(Rossow and Garder 1993) differs from the HBTM in 
several respects. It uses 3-hourly, 8 -km data sampled 
every 32 km and has slightly different criteria for 
selecting clear scenes. However, both techniques should 
yield similar clear-sky reflectances. Reflectances are 
used instead of albedos because the ISCCP assumes 
Lambertian reflectance for all land surfaces. The 
monthly mean hourly clear-sky reflectances are plotted 
in figure A3 for the four 2.5° regions near the comers of 
the SGP large-scale domain. The 1994 values closely 
track the 1985 HBTM results except during midmorning. 
The 1994 albedos in the upper right comer are slightly 
greater than those observed in 1985. The ISCCP results 
are consistent with the 1985 HBTM values except for the 
lower right quadrant. The ISCCP values are probably in 
error for this region because the morning values in the 
other regions are less than the afternoon values. This 
diurnal pattern is unlikely to change for the lower right 
region. 

Overall, the mean difference between the hourly 
1985 HBTM and 1994 ARM reflectances is -0.006 with 
an rms difference of +0.013 or 9 percent. The rms differ- 
ence between the regional means (integrated over the 
diurnal cycle), 0.0036 or 2.5 percent, is well within the 
ERBE interannual differences. If the dashed lines in fig- 
ure Al were used instead of equation (A5), the mean dif- 
ference would have been 7 percent, or almost twice the 
mean value obtained with equation (A5). The rms differ- 
ences between the hourly 1985 HBTM and ISCCP albe- 
dos are 8.6 percent or 7.6 percent if the lower right 
quadrant is excluded. Similarly, the 1994 LBTM and 
1985 ISCCP rms albedo differences are 11.1 percent and 
7.9 percent if the lower right region is excluded. The 
monthly rms differences between the daily mean 1994 
and ISCCP albedos is 7.8 percent. The same quantity 
computed for the ISCCP and HBTM albedos is 6.9 per- 
cent. From the consistency between the 1985 HBTM and 
1994 LBTM clear-sky albedos and for most ISCCP data, 
it is concluded that the GOES calibration procedure used 
for the 1994 data will produce reasonably accurate VIS 
albedos. 
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Figure Al. Intercalibration of April 1994 NOAA-11 AVHRR VIS radiances corrected to GOES angles and squared 
GOES-7 VIS counts; original regression fit shown by small dashed lines; final fit only used points for D 2 < 10000. 
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Figure A2. Mean April 1985-1989 and annual ERBE SW clear-sky albedos for 2.5° regions between 32.5°N and 
42.5°N and between 92.5°W and 105°W. 
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Appendix B 

Format of Archived ARM Satellite-Derived 
Cloud Properties 

The derived cloud properties are archived at NASA 
Langley Research Center. The Atmospheric Radiation 
Measurement (ARM) satellite-derived cloud dataset can 
be accessed on the Internet via the World Wide Web at 
the following uniform resource locator: 

http://albedo.larc.nasa.gov: 1123/arm. html 

The archived data are formatted in Network Common 
Data Form (netCDF) with version 2.3.2 of the netCDF 
software obtained from the Unidata Program Center. 
NetCDF interface software includes system callable utili- 


ties and input/output functions callable from C or 
FORTRAN. The system utility ' ncdump' can be used 
to generate a netCDF Common Data form Language 
(CDL) file that describes the format of data stored in a 
netCDF data file with or without including the variable 
data contained in the data file. This appendix contains a 
CDL file generated with the following commands for 
ncdump: 

ncdump ba$e__time, time__off set, latitude , 
longitude , level , view 

sgpgoes 7minni sXl . cl . 9404 05 .133000 

This sample file is a netCDF data file, 
sgpgoes 7minnisXl . cl .940405.133000, con- 
taining cloud products from 5 April 1994 over the South- 
ern Great Plains ARM region: 


netcdf sgpgoes 7minnisXl .cl . 940405 . 133000 { 
dimensions : 

latitude -20 ; 
longitude =28 ; 
level = 4 ; 
view = 2 ; 

time = UNLIMITED ; // (20 currently) 

variables : 

long base_time ; 

base_time : string = ”0:00:00 GMT 05 April 1994 " ; 

base_time : long_name = "base time in epoch” ; 

base__time : units = "seconds since 0:00:00 GMT 01 January 1970" ; 
double time_off set (time) ; 

time__of f set : long_name = "time offset from base time” ; 
time_off set : units = "seconds” ; 
float latitude (latitude) ; 

latitude: valid_range = -90. f, 90. f ; 
latitude : long_name = "north latitude" ; 
latitude: units ~ "degrees" ; 
float longitude (longitude) ; 

longitude :valid_range = -180. f, 180. f ; 
longitude : long_name = "east longitude" ; 
longitude : units = "degrees" ; 
float level (level) ; 

level : valid_range = 1 . f , 4 . f ; 
level : long_name = "cloud level" ; 
level: units = "unitless" ; 


level rvalue 

1 = "A 

value 

of 

1 . 

corresponds 

to 

low 

clouds . " 

level : value 

_2 = ”A 

value 

of 

2. 

corresponds 

to 

mid 

clouds . 11 

level : value 

_3 = "A 

value 

of 

3. 

corresponds 

to 

high 

l clouds. 

level : value 

_4 = "A 

value 

of 

4. 

corresponds 

to 

all 

clouds . " 


float view (view) ; 

view : valid_range = 1 . f , 2 . f ; 

view: long_name = "scene description" ; 

view: units = "unit less" ; 


20 



view: value_l = "A value of 1, corresponds to a clear scene.” ; 
view: value__2 = "A value of 2. corresponds to the total scene." 
float Cloud_Amount (time, level, latitude, longitude) ; 

Cloud_Amount : valid_jrange - O.f, 100. f ; 

C 1 oud_Amoun t : _Fi.l lvalue = - 8 8 8 . f ; 

Cloud_Amount :missing__value = -999. f ; 

Cloud__Amount : long_name = "cloud amount" ; 

Cl oud_Amount : units = "percent cloudy" ; 
float Visible_Optical_Depth (time, level, latitude, longitude) ; 
VisiblejDptical_Depth: valid_range = O.f, 200. f ; 

Vi s ib 1 e_Opt i c a 1 JD ep t h : __F i 1 1 Value - -888. f ; 
Visible_Optical_Depth:missing_value = -999 . f ; 
Visible_Optical_Depth:long_name = "visible cloud optical depth 
Visible_Optical_Depth: units = "unit less" ; 
float IR_Optical_Depth (time, level, latitude, longitude) ; 
IR_Optical_Depth: valid_range = O.f, 100. f ; 
IR^Optical^Depth-.^FillValue = -888. f ; 
IR__Optical__Depth:missing_value = -999. f ; 

IR_Optical_Depth: long_name - "IR cloud optical depth" ; 
IR_Optical_Depth : units = "unitless" ; 
float Emissivity (time, level, latitude, longitude) ; 

Emissivity : valid^range = O.f, l.f ; 

Emissivity :_FillValue = -888. f ; 

Emissivity :missing_value = -999. f ; 

Emissivity : long__name = "IR beam cloud emissivity" ; 

Emissivity : units = "unitless" ; 
float Cloud__Center__Height (time, level, latitude, longitude) ; 
Cloud_Center_Height : valid_range = O.f, 20 . f ; 
Cloud_Center_Height :_FillValue - -888. f ; 

Cloud_Center_Height :missing_value = -999. f ; 
Cloud_Center_Height : long_name = "cloud center height" ; 
Cloud_Center_Height : units = "kilometers" ; 
float Cloud__Top_Height (time, level, latitude, longitude) ; 
Cloud__Top_Height : valid_jrange = 0 . f , 20 . f ; 

Cloud_Top_Height :_FillValue = -888. f ; 

Cloud__Top__Height :missing_value = -999 . f ; 

Cloud_Top_Jieight : long_name = "cloud top height" ; 

C 1 oud__T op_He i ght : units = "kilometers" ; 
float Cloud__Temperature (time, level, latitude, longitude) ; 
Cloud__Temperature : valid__range = 160. f, 330. f ; 
Cloud_Temperature:_FillValue = -888. f ; 

C 1 oud_Tempe r at ur e : mi s s i ng_value = -999. f ; 

Cloud_Temperature : long_name = "uncorrected cloud temperature" 
Cloud_Temperature: units = "Kelvin" ; 
float Cloud__Thickness (time, level, latitude, longitude) ; 
Cloud__Thickness : valid_range - O.f, 20 . f ; 

Cloud__Thickness :__FillValue - -888. f ; 

Cloud__Thickness :missing__value = -999. f ; 

Cloud__Thickness : long_name = "cloud thickness estimate" ; 
Cloud__Thickness : units = "kilometers" ; 
float Reflectance (time, level, latitude, longitude) ; 

Reflectance :valid_range = O.f, 1 . 5f ; 

Reflectance :_FillValue - -888. f ; 

Ref lectance : mi s s 1 ng__value = - 9 9 9 . f ; 

Reflectance : long__name = "cloud narrowband reflectance at TOA" 



"unitless" 


Reflectance runits = 
float Albedo (time, level, latitude, longitude) ; 

Albedo :valid_range = O.f, l.f ; 

Albedo r^FillValue = -888. f ; 

Albedo rmissing_value - -999. f ; 

Albedo : long_name = "cloud narrowband albedo" ; 

Albedo runits = "unitless" ; 

float Cloud__Center_Temperature (time, level, latitude, longitude) ; 
Cloud_Center__Temperature: valid_range = 160. f, 330. f ; 
Cloud_Center__Temperature :_FillValue = -888. f ; 
Cloud_Center_Temperature r missing_value = -999. f ; 

Cloud___Center_Temperature r long__name = "cloud center temperature" ; 
Cloud__Center_Temperature runits - "Kelvin" ; 
float Cloud__Top_Temperature (time, level, latitude, longitude) ; 
Cloud_Top___Temperature r valid_range = 160.f, 330. f ; 
Cloud_Top_Temperature r_FillValue = -888. f ; 

Cloud_Top_Temperature rmissing_value = -999. f ; 

Cloud_Top__Temperature : long_name = "cloud top temperature" ; 
Cloud_Top__Temperature runits = "Kelvin" ; 
float Visible_Optical__Depthj3D (time, level, latitude, longitude) ; 
Visible__Optical_Depth__SD: valid_range = O.f, 100. f ; 
Visible_Optical_Depth_SD :_FillValue = -888. f ; 
Visible_Optical_Depth_SD rmissing___value = -999. f ; 

Visible_Optical_Depth_SD r long_name = "visible optical depth standard 
deviation" ; 

Visible_Optical_DepthJ3D runits = "unitless" ; 
float Cloud_Center__Temperature_SD (time, level, latitude, longitude) ; 
Cloud_Center_Temperature_SD r valid_range = O.f, 20 O.f ; 
Cloud_Center_Temperature__SD r_FillValue = -888. f ; 
Cloud__Center_Temperature_j3D rmissing_value = -999. f ; 
Cloud__CenterJTemperature__SD : long_name - "cloud center temperature 
standard deviation" ; 

Cloud__Center_Temperature___SD runits = "Kelvin" ; 
float Broadband_LW__Flux (time, view, latitude, longitude) ; 
Broadband_LW_Flux : valid_range = O.f, 400. f ; 
Broadband_LW_Fluxr_FillValue = -888. f ; 

Broadband_LW_Flux :missi.ng__value = -999. f 

Broadband_LW_Flux: long_name = "broadband LW flux for clear and total 

scene" ; 

Broadband_LW__Flux runits - "watts per square meter" ; 
float Narrowband_IR__Flux (time, view, latitude, longitude) ; 
Narrowband_IR___Fluxr valid_range = O.f, 10 O.f ; 

Narrowband_IR_Flux r_FillValue = -888. f ; 
Narrowband_IR_Flux:missing_value = -999. f ; 

Narrowband_IR_Fluxr long_name •— "narrowband IR flux for clear and 
total scene" ; 

Narrowband__IR_Flux : units = "watts per square meter" ; 
float Broadband__SW_Albedo (time, view, latitude, longitude) ; 
Broadband_SW__Albedo : valid_range = O.f, l.f ; 

Broadband__SW_Albedo :_Fil lvalue = -888. f ; 
Broadband__SW_Albedormissing__value = -999. f ; 

Broadband__SW_Albedo r long_name = "broadband SW albedo for clear and 
total scene" ; 

Broadband SW Albedo runits = "unitless” ; 
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float Narrowband_VIS_Albedo (time, view, latitude, longitude) ; 
Narrowband__VIS_Albedo : valid_range - O.f, l.f ; 
Narrowband_VIS_Albedo :__Fi 11 Value = -888 . f / 

Nar rowband_VI S__Albedo : mi s s ing_value = -999. f ; 

Narrowband_VIS_Albedo : long_name = "narrowband VIS albedo for clear 
and total scene” ; 

Narrowband_VIS_Albedo : units = "unitless” ; 
float ClearJTemperature (time, latitude, longitude) ; 

Clear_Temperature : valid__range = 160. f , 330 . f ; 

Clear_Temperature :_FillValue = -888. f ; 

Clear_Temperature :missing_value = -999. f ; 

ClearJTemperature : long_name = "clear IR temperature" / 
ClearJTemperature : units = "Kelvin” ; 
float Clear jTemperature_SD (time, latitude, longitude) ; 
Clear_Temperature_SD : valid_range ■= O.f, 20 O.f ; 
Clear_TemperaturejSD : jFillValue = -888. f ; 
Clear_TemperaturejSD:missing_value = -999. f ; 

Clear_Temperature_SD : long_name = "clear IR temperature standard 
deviation" ; 

Clear_TemperaturejSD : units = "Kelvin" ; 
float NarrowbandjVIS_AlbedOjSD (time, latitude, longitude) ; 

Nar r owbandjVI S_Albedo_SD : valid_range - 0 . f , l.f ; 
Narrowband__VIS_Albedo_SD :_FillValtie = -888. f ; 
Narrowband_VIS_Albedo_SD :missing_value = -999. f ; 

Narrowband_VIS_Albedo_SD : long__name = "clear narrowband VIS albedo 
standard deviation" ; 

Narrowband_yiS_Albedo_SD : units = "unit les s " ; 
float Clear__VIS_Reflectance (time, latitude, longitude) ; 

Clear_VIS__Ref lectance : valid_range -O.f, l.f; 

Clear_VIS__Ref lectance : jFillValue = -888. f ; 

Clear_VISjRef lectance :missing_value = -999. f ; 

Clear jVIS_Ref lectance : long_name = "clear narrowband VIS reflectance" 
Clear_JVIS_Ref lectance : units = "unitless" ; 
float Average jTotal_Temperature (time, latitude, longitude) ; 
Average__Total_Temperature : valid_range = 160. f, 330 . f ; 
Average_Total_Temperature : ^FillValue = -888. f ; 
Average_Total_Temperature :missing_value = -999. f ; 
Average_Total_Temperature : long_name = "total scene average 
temperature" ; 

Average_Total_Temperature : units = "Kelvin" ; 
float SolarjZenith__Angle (time, latitude, longitude) ; 

Solar_ZenithjAngle: valid^range = O.f, 90. f ; 

Solar_ZenithjAngle: jFillValue = -888 . f ; 

Solar_Zenit hjAngle : mi s s ing_value - -999. f ; 

Solar_Zenith_Angle : long_name = "solar zenith angle" ; 
Solar_Zenith_Angle : units = "degrees" ; 
float Viewing_Zenith_Angle (time, latitude, longitude) / 
Viewing_Zenith_Angle : valid_range = O.f, 90. f ; 
Viewing_Zenith_Angle: jFillValue = -888. f ; 

Viewing_Zenith_Angle zmissing^value = -999 . f ; 

Viewing jZenith_Angle : long_name = "viewing zenith angle" ; 
Viewing_Zenith_Angle : units = "degrees " ; 
float Relative_Azimuth_Angle (time<, latitude, longitude) ; 

Re lative__Az imuth__Angle : vali derange = O.f, 180. f ; 

Relat ivejAzimuth__Angle : _FillValue = -888 . f ; 



Re 1 at i ve_A z imu t h_An gle:missin g__va 1 u e = -999 .f ; 

Relative_Azimuth_Angle : long_name = "relative azimuth angle” ; 

Re 1 at i ve_Az imut h_Angl e : units = "degrees” ; 
float alt ; 

alt : long__name = "Dummy altitude for Zeb" ; 
alt : units = "unitless" ; 
float lat ; 

lat : long_name = "Northernmost north latitude for Zeb" ; 
lat: units = "degrees" ; 
float Ion ; 

Ion : long_name = "Westernmost west longitude for Zeb" ; 

Ion: units = "degrees" ; 

// global attributes: 

: title = "0.5 degree LBTM cloud products derived from GOES for ARM 
Great Plains" ; 

: source = "NASA Langley Research Center" ; 

: version = "LBTM ARM 1.0.0" ; 

: netCDF = "net CDF 2.3.2" ; 

: reference = "Minnis, P., Heck, P . W. , and Young, D. F . , 1993: Infer- 
ence of Cirrus Cloud Properties Using Satellite-observed Visible and Infrared 
Radiances. Part II: Verification of Theoretical Cirrus Radiative Properties. J. 
Atmos . Sci . , 50, 1305-1322 . " ; 

: visible_calibration = "The visible radiance was calculated according 
to R ± a * D * D + b, where R is the visible radiance, D is eight bit counts, 
a = 0.0127, and b=-.28." ; 

: inf rared_calibration = "The nominal GOES calibration was used for 

infrared." ; 

: shortwave_NB/BB_correlation ■■= "The shortwave narrowband/broadband 
correlation is given by Ab = Abclr * ( 1 - X ) + Abe Id * X, with Abclr = a + b * 
Anclr + c * In ( 1 / uo ) and Abcld - d + e * Ancld + f * Ancld * Ancld + g * In ( 

1 / uo ), where Ab is broadband albedo, Abclr is clear sky broadband albedo, 
Abcld is cloudy sky broadband albedo, X is cloudy sky scene fraction, Anclr is 
clear sky narrowband albedo, Ancld is cloudy sky narrowband albedo, uo is the 
cosine of the solar zenith angle, a = 0.1218, b = 0.3842, c = 0.0605, d = 0 . 0588, 
e = 0.8623, f = -.1190, and g = 0.0624." ; 

: shortwave_NB/BB_reference = "Doelling, D . R. , Young, D. F . , Arduini, 
R. F . , Minnis, P . , Harrison, E. F . , and Suttles, J. T. , 1990: On the Role of Sat- 
ellite-measured Narrowband Radiances for Computing the EarthVs Radiation Bal- 
ance. Proc. Seventh Conference on Atmospheric Radiation, San Francisco, CA, July, 
155-160." ; 

: longwave_NB/BB_correlation = "The longwave narrowband/broadband cor- 
relation is given by Mb = a + b * Mn + c * Mn * Mn + d * Mn * In ( h ) , where Mb 
is broadband flux, Mn is narrowband flux, h is the average relative humidity, in 
percent, above the GOES level, a = 63.6, b = 6.628, c = -.0278, and d = -.332." ; 

: longwave_NB/BB_reference = "Minnis, P . , Young, D. F . , and Harrison, 
E. F . , 1991: Examination of the Relationship between Outgoing Infrared Window and 
Total Longwave Fluxes Using Satellite Data. J. Climate, 4, 1114-1133." ; 

data: 

base_time = 765504000 ; 

time_of f set = 48600, 50400, 52200, 54000, 55800, 57600, 61200, 63000, 64800, 
66600, 68400, 70200, 72000, 73800, 75600, 77400, 79200, 81000, 82800, 

84600 ; 
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latitude = 41.75 , 41.25 , 40.75 , 40.25 , 39.75 , 39.25 , 38.75 , 38.25 , 
37.75 , 37.25 , 36.75 , 36.25 , 35.75 , 35.25 , 34.75 , 34.25 , 33.75 , 
33.25 , 32.75 , 32.25 ; 

longitude = -104.75 , -104.25 , -103.75 , -103.25 , -102.75 , -102.25 , 
-101.75 , -101.25 , -100.75 , -100.25 , -99.75 , -99.25 , -98.75 , 
-98.25 , -97.75 , -97.25 , -96.75 , -96.25 , -95.75 , -95.25 , -94.75 , 
-94.25 , -93.75 , -93,25 , -92.75 , -92.25 , -91.75 , -91.25 ; 

level =1,2*3, 4 ; 

view =1,2 ; } 
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Figure 1 . Satellite data analysis grids and surface elevation map for ARM Southern Great Plains locale. 
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Figure 2. Mean hourly surface shelter and GOES-derived clear-sky temperatures for April 1994 over ARM SGP meso- 
scale domain. 
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Figure 3. Observed and regression fit of A7' Vf;i for April 1994 over ARM SGP mesoscale domain. 
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(a) Clear-sky data. 
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(b) AH data. 

Figure 5. Correlations of ERBS SW and GOES VIS albedos for April 1985; curves indicate regression fits at particular 
solar zenith angles. 
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Figure 6. Correlations of ERBS LW and GOES IR fluxes for April 1985; curves indicate regression fits at particular rel- 
ative humidities. 
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Figure 7. Cloud optical depths overlaid on GOES VIS image for 1930 UTC, April 14, 1994. 



Figure 8. Total cloud amounts overlaid on GOES IR image for 1930 UTC, April 14, 1994. 
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Figure 9. GOES VIS and IR images overlaid with ARM SGP mesoscale domain for 1800 UTC, April 25, 1994. 
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Figure 10. Cloud amounts for 1800 UTC, April 25, 1994. 
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Figure 11. Cloud optical depths for 1800 UTC, April 25, 1994. 
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Figure 12. High and total cloud center and top heights for 1800 UTC, April 25, 1994. 
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Figure 13. Cloud thickness for 1800 UTC, April 25, 1994. 
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Figure 14. VIS and SW albedos for 1800 UTC, April 25, 1994. 
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Figure 15. IR temperatures and LW fluxes for 1800 UTC, April 25, 1994. 
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Figure 18. Mean half-hourly cloud amounts, optical depths, and heights for 0.3°, 0.5°, and 2.5° regions centered on SCF 
for April 5 to May 1, 1994. 
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Figure 19. Frequency of cloud-top height occurrences for 0.5° region (36.75°N, 97.25°W) that includes SGP Central 
Facility. 
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Figure 20. Layer cloud thickness for 0.5° region (36.75°N, 97.25°W) that included SGP Central Facility. 



0 1 2 3 4 5 6 7 

CLOUD THICKNESS (km) 



45 



42° N 


Total Cloud Amount 




105° W 91 ° W 

Mid Cloud Amount 


32° N 



42° N 


Low JC loud Amount 


32° N 
105° W 



91° W 


105° W 


91° W 



0 10 20 30 40 50 60 70 80 

Cloud Amount (%) 

Figure 21. Mean cloud amount for 1330 to 2330 UTC for April 5 to May 1, 1994. 
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Figure 22. Mean cloud optical depths for 1330 to 2330 UTC for April 5 to May 1, 1994. 
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Figure 23. Mean cloud-center heights for 1330 to 2330 UTC for April 5 to May 1, 1994. 
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Figure 24. Mean cloud thickness for 1330 to 2330 UTC for April 5 to May 1, 1994. 
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Figure 25. Mean VIS and SW albedos for 0.3° and 2.5° (average of 0.5°) regions centered on SCF for April 5 to May 1, 
1994. 
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Figure 27. Mean VIS and SW albedos for 1330 and 2330 UTC for April 5 to May 1, 1994. 
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Figure 28. Mean IR temperatures and LW fluxes for 1330 to 2330 UTC for April 5 to May 1, 1994. 
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(a) 0.5° centered on WAFB at 38.75°N, 93.75°W; 1.0° centered on WAFB at 38.73°N, 93.55°W. 
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(b) 1.0° centered on TAFB at 35.42°N, 97.38°W. 

Figure 29. Mean LBTM cloud amounts for each tenth of ground observed sky cover for April 5 to May 1, 1994; missing 
values filled by using linear interpolation for the satellite analysis regions. 
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(b) 1 .0° centered on TAFB at 35.42°N, 97.38°W. 

Figure 3 1 . Mean hourly matched cloud amounts from surface and from satellite for April 5 to May 1 , 1994. 
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